Observing fermionic statistics with photons in arbitrary processes 
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Quantum mechanics defines two classes of particles — bosons and fermions — whose exchange statis- 
tics fundamentally dictate the dynamics of quantum systems. Here we develop a scheme that uses 
multi-partite entanglement to directly observe the correlated detection statistics of any number of 
fermions in any physical process. This approach relies on sending each of the entangled particles 
through identical copies of the process and by controlling a single phase parameter in the entangled 
state, the correlated detection statistics can be continuously tuned between bosonic and fermionic 
statistics. We implement this scheme via two entangled photons shared across the polarisation 
modes of a single photonic chip to directly mimic the fermion, boson and intermediate behaviour of 
two-particles undergoing a continuous time quantum walk with an average similarity of 93.6±0.2% 
with ideal models. The ability to simulate fermions with photons is likely to have applications for 
verifying computationally hard boson scattering in linear optics and for observing particle correla- 
tions (including fractional statistics) in analogue simulation using any physical platform that can 
prepare the entangled state prescribed here. 



A particularly striking departure from classical physics 
in quantum mechanics is the existence of two fundamen- 
tal particle classes characterised by their behaviour under 
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pairwise exchange (aja[ = ±0^0]). For example an arbi- 
trary number of bosons (obeying Bose-Einstein statistics: 
'-I-') can occupy a single quantum state, leading to bunch- 
ing [T] and Bose-Einstein condensation; whereas fermions 
(obeying Fermi-Dirac statistics: '— ') are strictly forbid- 
den to occupy the same state via the Pauli-exclusion prin- 
ciple Behaviour intermediate between bosons and 
fermions can also be considered; particles confined to two 
dimensions known as anyons [3j are predicted to exhibit 
intermediate behaviour with fractional exchange statis- 



tics (ai-ajj, = e^'^'aj.aj, Q < 4> < t) — a concept applied in 
superconductivity [3] and used to explain the fractional 
quantum Hall effect [S]. 

Here we show that by tuning one phase parameter, 
entanglement and a single experimental setup can ex- 
actly simulate the measurement statistics of an arbi- 
trary mode transformation A (Fig. [T](a)) of N identical 
fermions, bosons or particles with fractional statistics. 
It is well known that the symmetric and anti-symmetric 
states IV'^) = \H)^ \ V)2^\V)^ \H)2 of two polarised pho- 
tons input into a 50% reflectivity beam splitter exhibits 
fermion-like anti-bunching and boson-like bunching, re- 
spectively 9J. Furthermore, intermediate behaviour for a 
beamsplitter has been observed by controlling a phase to 
tune between the two Bell states \^^) and |^~) PUIITT] . 
It is less widely appreciated that the resulting detection 
statistics can also be exactly mimicked with two pho- 
ton path- number entangled NOON states [T^]; however 
this behaviour does not generalise to processes other than 
beamsplitters, the simplest example being the identity 
operation for which bunching always occurs regardless of 
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FIG. 1: Simulating quantum interference of quantum particles, 
(a) Quantum interference of identical particles with arbitrary 
statistics launched into a network described by matrix A leads to 
quantum correlated detection events at the output, (b) Experi- 
mental simulation of quantum interference can be directly achieved 
by sharing Af-partite, A''-level entanglement across N copies of A, 
where the choice of simulated particle statistics is controlled by the 
phase of initial entanglement. 



the phase of the NOON state. The approach we present 
shares A^-partite, A'^-level entanglement, parameterised 
by phase 0, across N copies of A (Fig. [l] (b)). 
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After presenting the scheme in the case of two parti- 
cles in an arbitrary linear evolution, we experimentally 
demonstrate the simulation of arbitrary particle statis- 
tics with a complex photonic network — a continuous time 
quantum walk [51 [7] .We achieve the required entangled 
state with post-selected polarisation entanglement shared 
across the TE and TM modes of an evanescently cou- 
pled waveguide array. This enables us to directly observe 
the evolution of symmetric and anti-symmetric two par- 
ticle wave functions in a coupled oscillator Hamiltonian 
by switching the phase of the entangled state between 
and TT. Mixed-symmetric quantum walk behaviour is 
also observed by tuning the entangled state's phase be- 
tween and TT, corresponding to behaviour intermedi- 
ate between bosons and fermions. We note the fermion 
simulation directly maps to the dynamics of two excita- 
tions in a spin— i chain consisting of ten sites: Using the 
Jordan- Wigner transformation, excitations correspond to 
fermion creation operators and the vacuum to the ground 
state of each spin [8j . We finally provide a mathematical 
demonstration of our scheme in the case of an arbitrary 
number of particles evolving through any linear Hamilto- 
nian. It is important to note that neither the two-particle 
case, nor the N-particle generalization, require any post- 
selection and are, in principle, scalable. 



Results 



detected at outputs r and q) is given by 



(4) 



for bosons [14. {((> — 0) and fermions [13j {(p — tt). Here, 
the process A is represented by a matrix whose element 
Arj is the complex transition amplitude from j to r. 
Non-classical behaviour is evident in the constructive and 
destructive interference between typically complex terms 
inside the |-|^, dictated by = uj9. By defining fractional 
exchange algebra (l][3) for < 6* < tt and some integer 
winding number w, we can observe the effect of fractional 
exchange statistics on quantum interference between two 
indistinguishable anyons undergoing braiding operations 
dependent upon A; Eq. Q represents the resulting quan- 
tum interference between two abelian anyons (0 < < tt) 
in superposition of braiding tu times and not braiding at 
all. To see this we define without loss of generality the 
input labels hy j < k and the output labels hy r < q and 
consider two identical anyons created at modes j and 
fc, transformed according to A and detected at modes r 
and q. Correlated detection at r and q reveals two in- 
distinguishable possibilities that correspond to the two 
terms in Eq. Q: A maps the anyon at j to mode r and 
the anyon at k to mode q with amplitude Aj- jAq ^] or 
A maps the anyon at j to mode q and the anyon at k 
to mode r via uj clockwise braids with fractional phase 
and amplitude AqjAr^k- From this definition, we use 



relations (l][3) to derive expression ^ for < 9 < TT and 
(j) — Lo9 (see Supplementary Information). 



Exchange statistics refer to a wave-function acquiring 
phase after interchanging two indistinguishable quan- 
tum particles. Occupation of labeled modes j are mod- 
elled by annihilation operators aj and creation operators 

that obey the commutation relations 

«]4-e'M«l = (1) 

ajttk — e'''^akaj = (2) 
«i4 - e"^4aj = Sj,k (3) 

for j < k. In three dimensions particles are either bosons 
{(j) — 0) or fermions {(f) = tt) and successive swapping 
is represented by the permutation group while the cor- 
related detection statistics are associated respectively to 
the permanent and determinant of a given matrix |13j . 
In one- and two-dimensions anyons [5] can exist where 
swapping two identical anyons to times (represented by 
the braid group) results in the combined wave-function 
acquiring phase (f) — uj9 where < < tt and the positive 
(or negative) winding number lo keeps track of the num- 
ber of times two anyons braid in trajectory space with 
clockwise (or anti-clockwise) orientation 

The quantum correlation function for two identical 
non-interacting particles (input into j and k of A, and 



Simulation of two-particle statistics 

To physically simulate quantum interference given by 
Eq. Q, we use two copies of A, denoted A°- and A^, to- 
gether with bi-partite entanglement of the form |'0(0)) — 
^'^'^ °^k^\) |0) 1^ where creation operators and 
h\ are associated to identical networks A" and A^ respec- 
tively; creation operator sub-indices denote the modes on 
which the creation operators act. Once |V'(0)) is gener- 
ated, the class of input entangled particles becomes irrel- 
evant; only one particle is present in each of A" and A^ , 
therefore no further quantum interference between parti- 
cles will occur. This implies application with other quan- 
tum platforms, using for example entangled fermions to 
simulate Bose-Einstein quantum interference. Evolution 
of the two particle state through the tensor product of 
transformations is therefore modelled by 

\m) E ( ^°^-^^'^-'"^'^°"'''^'^ ) albl 10) (5) 

It follows that correlated detection probability across 
outputs r of A°- and q of A^ simulates Eq. Q accord- 
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ing to 

s,t 

1 n 1 

— -14*^ A- p^"^ A'^ a'' r = -T*^ (fW 

— 2 \ "^-i ^r-,fc^qj| — 2 t-,? V"J 

where the 1/2 pre- factor is due to normahsation of the 
initial entanglement. Since the phase of the entangled 
state can be defined as 4> = uj9 , this method simulates the 
interference statistics of two anyons of fractional phase 9 
that have not braided, with two anyons also characterised 
by 6 that have braided lj times (equivalently anyons that 
are in a superposition of braiding loi and uj2 times for 
(j> = (cji — 102)0 — the ideal statistics are the same). 

Experiment 

We experimentally simulated Eq. Q for arbitrary (j) 
in a coupled oscillator Hamiltonian equivalent to sharing 
entanglement across two continuous time quantum walks, 
realised with arrayed waveguides [5] (see Supplementary 
Information). We measured the correlations of the po- 
larisation entangled state {\H). \V)f^ + e'^ \V)j \H)^)/^/2 
shared across two copies of a process A: A° and A'^ , re- 
alised as the TE and TM modes of a waveguide array 
fabricated in silicon oxynitride (SiO^^Ny) (Fig. [2](b), see 
also Supplementary Information). The waveguides and 
butt-coupled polarisation maintaining fibres (PMFs) are 
birefringent to ensure negligible cross-talk between hori- 
zontal and vertical polarisations. Detection statistics of 
the output photons in the , basis were anal- 
ysed using polarising beam splitters (PBS) and pairs of 
multi-mode fibre coupled avalanche photo-diode single 
photon counting modules (SPCM) (Fig. [2] (c)). By in- 
putting the state {\H)^ + \V).){\H)^^ + e"^ \V)^)/2 and 
detecting one \H) photon and one \V) photon at any 
pair of outputs of A° and A'^ we post-selected the input 
state {\H)^ -he*('^+'=) \V) ^ \H)^)/y/2, with a measur- 
able phase off-set e (from birefringence of PMF and the 
waveguide circuit). This phase e was compensated us- 
ing a polarisation phase rotation (Fig. [2] (a)). Note that 
because of this compensation, the approach described 
is immune to birefringence of the 2m PMF and input 
waveguides. 

We measured two sets of correlation matrices for two 
different inputs: immediate neighbouring waveguides in 
the centre of the array (j, fc) — (—1,0) (Fig.jsj); and next 
but one neighbouring waveguides (j,^ = (~li 1)j with 
vacuum input in waveguide (Fig. W|. Each en- 
try is the probability to detect a horizontal and vertical 
polarised photon at outputs r and q, respectively. Five 
cases are shown for = 0, 7r/4, 7r/2, 37r/4, tt (left: (a), (c), 
(e), (g), (i)). The corresponding ideal correlation matri- 
ces given by Tf ,^ in Eq. m are also plotted (right: (b). 




FIG. 2; Experimental setup, (a) The parametric down-conversion 
based photon pair source with quarter- waveplates (QWP), half- 
waveplates (HWP) and a polarisation beamsplitter cube (PBS) 
used as a polariser, to control the initial state before input into the 
chip. Down conversion in the nonlinear bismuth borate (BiBO) 
crystal yields degenerate photon pairs (808nm, filtered with 2nm 
full width, half maximum interference filters, IF) collected into po- 
larisation maintaining fibres (PMF). See Supplementary Informa- 
tion for further details, (b) Two copies of the quantum walk unitary 
are realised by accessing the TE and TM modes of the polarisa- 
tion preserving waveguide array. Three input modes of the array 
are used in the experiment (-1,0,1) accessed via waveguide bends 
that fan to a pitch of 250/^m and butt-coupled to an array of PMF. 
Output is accessed via a fan of waveguide from the coupling region 
to a pitch of 125/im. See Supplementary Information for details. 

(c) The detection scheme consists of five PBS and single photon 
counting module (SPCM) pairs. PMF coupled to the chip launch 
the output photons onto each PBS which separates TE and TM 
modes guided in the waveguide. The two outputs of each PBS 
are collected into multi-mode fibre (MMF) coupled to SPCMs for 
correlated detection via counting electronics. 

(d) , (f), (h), (j)). For the {j,k) - (-1,0) input (Fig. § 
we observe the expected bunching behaviour of bosons 
for (j) — Q which is contrasted with the observation of the 
Pauli-exclusion in the diagonal elements of the matrices 
for (f> — TT. We note that a form of exclusion on the diag- 
onal entries would also occur for experiments with 
hard core bosons, but the remaining correlation matrix 
would still be dependent upon Bose-Einstein exchange 
statistics. The spectrum of phase < < tt exhibits a 
variable Pauli-exclusion and a continuous transition be- 
tween simulation of Bose-Einstein and Fermi-Dirac quan- 
tum interference. In the bosonic case, we observe a less 
pronounced bunching behaviour for the {j,k) = (—1,1) 
input (Fig. [4] (a)) than for the (j, fc) = (—1,0) input 
(Fig. [3] (a)), as expected [T^. This behaviour is inverted 
in the fermionic case, since the particels tend to propa- 
gate to the same peak of the walk (while still not occu- 
pying the same waveguide) for the {j, k) = (—1, 1) input 
(Fig. [4] (i)), and tend to propagate in different peaks for 
the (j, fc) = (-1,0) input (Fig. [3] (i)). This is the result 
of the different phases in the Hamiltonian that describes 
the quantum walk evolution. 

Agreement between the ideal implementation of our 
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scheme, with perfect reahsations of A"" and A^, and 
the data in Figs. 3j4 is quantified in Table |l] using the 
similarity 5* between the ideal correlation matrix 
and the meas ured dis tribution of detection statistics P'^: 

S = (E.,, \/^t^,)VEr,,^t..,Er,,Pr^,,^ ^ geucralisa- 
tion of the average (classical) fidelity between two proba- 
bility distributions that ranges from 0% (anti-correlated 
distributions) to 100% for complete agreement. We also 
calculated S when we deliberately made the photons dis- 
tinguishable by temporally delaying one with respect to 
the other. This enabled physical simulation of the clas- 
sical correlations F*^ = -I- that 
are independent of particle statistics. For these mea- 
surements, we found agreement with the ideal model of 
S = 97 A ± 0.7% and S = 97.5 ± 0.6% for inputs (-1, 0) 
and (—1,1) respectively. These data indicate that A" 
and A^ are near identical while keeping crosstalk between 
horizontal and vertical polarisation to a minimum. While 
non-unit similarity {S — 100%) in Table |l] is attributed 
to non-perfect entanglement generation, the high level of 
agreement between the ideal and measured cases indi- 
cates accurate quantum simulation of arbitrary quantum 
interference in a uniform coupled oscillator Hamiltonian. 



Generalization to N particles 

Using multi-partite entanglement allows simulation of 
quantum interference of an arbitrary number N of identi- 
cal non-interacting particles. N indistinguishable parti- 
cles with exchange statistics (j), subjected to mode trans- 
formation A at inputs v = {z^i, 1/2, j^at}, are corre- 
lated due to quantum interference across outputs fl = 
{/^i,/^2, •■•,A*Af} according to 
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where Gi, is an element of the permutation group Sm act- 
ing on v and r {ui,) is the minimum number of neighbour- 
ing swaps that maps ay to v] a^j denotes the element j 
of Uy. For ^ = and (j) — n, Eq. ([t]) is the modulus- 
square of the permanent and determinant respectively of 
an X TV subset of A. For < (/)< tt Eq. ^ — which is 
the modulus-square of the imminent of the N x N subset 
of A — is derived in the Supplementary Information. The 
multi-particle correlation function Eq. Q is simulated by 
subjecting the generalised TV-partite, A^-level entangled 
state 
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(8) 



to A'' copies of A such that one particle a(j")^ is launched 
into each copy A^^\ For example, for inputs j,k,l and 
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FIG. 3: Experimental simulation of particle statistics correlation 
functions for input (—1, 0). Experimental correlation matrices (left: 
a,c,e,g,i) for simulating quantum interference of non-interacting 
particle pairs with five different particle statistics. Ideal correla- 
tion matrices for this experiment are plotted (right: b,d,f,h,j) us- 
ing Eq. (j4]l. The data ranges from corresponding to Bose-Einstein 
statistics {(p = 0, (a,b)) through fractional statistics {<j) = 7r/4, 
(c,d); = tt/2, (e,f); ij> = 37r/4, (g,h)) to Fermi-Dirac statis- 
tics (</) = TT, (i,j)). (k) The correlation matrices are displayed 
as normalised probability distributions with white for probability 
Pr',q = and black for P^ q = 1. Solid coloured squares represent 
measured correlations, hatched squares correspond to coincidences 
that cannot be measured due to mis-match between fibre array 
separation and waveguide pitch (see Supplementary Information). 
Coincidence rates are corrected for relative detection efHcicncy. 
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FIG. 4: Experimental simulation of particle statistics correlation 
functions for input (—1, 1) with data laid out as for Fig. [s] 



Phase, (j> 


S (%) for input (-1,0) 


S (%) for input (-1,1) 





95.8 ±0.6 


90.9 ±0.5 


7r/4 


95.2 ±0.7 


93.6 ±0.6 


tt/2 


94.1 ±0.7 


93.5 ±0.6 


37r/4 


94.4 ±0.7 


94.1 ±0.6 




91.6 ±0.7 


92.7 ±0.6 



TABLE I: Similarity S for each experimental simulation of particle 
quantum interference with parameter (p, for input (—1,0) given in 
Fig. [3] and input (—1, 1) given in Fig.|4] 



A'' = 3, the required entangled input state is equivalent 
to three entangled qutrits, shared across three copies of 
A denoted A^'^A^ and A": 



v/6( 



W'^alh\c] + e^'^alb]cl + e^'^ajblc]) |0) (9) 

iV-fold coincidence detection across outputs /ij of A^^^ 
equals -^^^p- Note as for the case with two particles, 
the 1/iV! pre- factor can be canceled by considering all 
A'^! permutations of the outputs, therefore directly sim- 
ulating quantum interference of N bosons (^ = 0) or N 
fermions {(j> — tt). For example, in the 2-particle case, the 
1 /2 pre-factor in Eq. ([6| is canceled by summing both the 
detection outputs of a particle at r on A", g on A'' and at 
r on A^, q on A"". All other values of (j) allow continuous 
transition between Fermi and Bose quantum interference. 
Note that the approach of initialising the simulation with 
a pre-symmetrized (or pre-antisymmetrized) state is in 
the spirit of first quantised quantum simulation [TS], for 
which quantum algorithms have been designed to effi- 
ciently generate the required states [El [17] (see Supple- 
mentary Information). In a photonic realization of the 
proposed sheme, the required N-particle entangled state 
can be generated using additional photons and feedfor- 
ward to implement quantum gates following the KLM 
scheme [18]. 



Discussion 

Correlations resulting from quantum interference of 
multiple indistinguishable non-interacting bosons in large 
linear optical circuits cannot be efficiently simulated with 
classical computers [T9U21j . Large scale repeated quan- 
tum interference is achievable with photons in waveguides 
[3 122H21] , which is likely to provide the platform for the 
realisation of such a boson scattering device, to realise 
an engineered system that can outperform the limits of 
classical computation. A major issue will be verifying 
boson scattering correlations that exceed the capabili- 
ties of conventional computers. In contrast, the anal- 
ogous situation for non-interacting fermions is efficient 
with classical computers |36j . The ability to tune con- 
tinuously between tractable and intractable phenomena 
with a single setup could be essential for verification of 
bosonic scattering. 

Other possible applications for simulating particle 
statistics with photonics include preparing appropriately 
symmetrized states in quantum chemistry simulations 
|15) and combinatorics, where exchange statistics have 
been shown to alter mathematical graph structure [2S] — 
three-dimensional waveguide architectures [26] may be 
important here. Our simulation of arbitrary fractional 
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exchange statistics also provides a novel method to study 
fractional quantum interference phenomena |27j ; whether 
our approach can be used in the context of topolog- 
ical fault tolerant quantum computation with photons 
[551 HH] is an open question. 

The method we report can be directly applied to any 
quantum process realisable with linear optics |30) and 
is not reliant on using polarisation — the entangled state 
could be encoded in time-bin or interferometrically stable 
path with reconfigurable waveguide circuits |31| . In order 
to introduce interactions between simulated particles, it 
may be possible to use effective non-linearities using aux- 
iliary photons [m [32] . In the long term, the method we 
have demonstrated here could be applied to other quan- 
tum platforms (consisting of fermions or bosons) where 
nonlinear behaviour may be more readily exploited. 
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Methods 

The photon pairs are generated in a 2mm thick, 
nonlinear bismuth borate BiBsOg (BiBO) crystal phase 
matched for non-colinear type-I spontaneous parametric 
downconversion, pumped by a 60 mW 404 nm continu- 
ous wave diode laser. The photons are spectrally filtered 
using 2nm full width at half maximum interference filters 
centred on 808nm and collected with aspheric lenses 
focusing on polarisation maintaining fibres (PMFs). 
The typical two-photon count rate directly from the 
source with two 60% efficient single photon counting 
modules (SPCMs) is 50kIIz. For a coherent input state, 
arrival time of the two photons at the input of the 
waveguide array is made identical via characterisation 
by performing a Hong-Ou-Mandel experiment |33]. The 
polarisation state \D) = {\H) + \V))/^ is then prepared 
for both photons using a half wave-plate with optic axis 
set to 7r/8 radians. The polarisation in the left arm is 
matched and purified to that of the right using a rotated 
polarisation beam splitter in transmission. The two 
photon state launched into the waveguide array is then 
{\H)j + \V)j)i\H)^ + e'^'^+'''> |F)fc)/2. The collecting PMF 
and the birefringence of the waveguide array preserve 
horizontal and vertical polarisation components for both 
photons and prohibit crosstalk between horizontal and 
vertical light; PMF-coupled polarisation beam splitters 



at the output of the waveguide array allow two-photon 
correlated detections across the two orthogonal po- 
larisations to project the initial input state to be the 
entangled state i\H)^ + e*W+^) \V)^ \H)^)/s/2- the 
extra phase e is compensated using a Z{0) rotation on 
the right arm of the photon source comprised of two 
quarter waveplates (optic axis at 7r/4 radians) with a 
half waveplate between (optic axis at — (tt + -I- e)/4 
radians, see Fig. 1 of the main text). 



The waveguide array was fabricated from silicon 
oxynitride (SiON) thin films were grown by plasma 
enhanced chemical vapour deposition (PECVD) onto 
thermally oxidised < 100 > oriented 100-mm silicon 
wafers applying an Oxford 80 Plasmalab Systems. A 
deposition process based on 2% Sill4/N2 and N2O was 
performed at 60W power, 187.5 kHz frequency, 650 
mTorr pressure and 300° C substrate temperature. The 
SiON films were annealed for 3hours at 1150°C in a 
nitrogen atmosphere. By optical inspection at 632 nm 
wavelength a film thickness oi d — 603 (±3) nm and a 
refractive index (for TE) of n^E = 1.5232 (±1 x lO'^) 
were measured. The channel waveguides were obtained 
from standard lithography followed by reactive ion 
etching in CHF3/O2 chemistry at 350 W plasma power. 
The etch depth [d = 613 ± 5nm) was confirmed by 
profilometry. A b-fiin thick upper cladding layer was 
deposited in two steps applying a Low Pressure CVD 
TEOS-process and standard PECVD Si02 deposition. 
The upper cladding was annealed at 1150°C for 3 hours. 
For a core width of 1.8-/xm the total birefringence of the 
waveguide, including material and form birefringence 
[35], is of the order An — titm — ^^te = 2 * lO""*. 
This value provides polarization preserving propagation 
for the TE and TM photons throughout the device. 
Moreover, the birefringence is small enough to ensure 
a nearly identical mode overlap, and coupling ratio, 
between different waveguides in the coupled array for 
TM and TE modes. This implies that the processes 
A°- = A^^ and — A^^'^ are matched, as required 
and verified in the reported experiment. The coupled 
waveguide arrays were accessed with lithographically 
patterned bends (minimum bend radius of 600/im) that 
gradually mapped the 250/Ltm pitched input (for 250/im 
pitched PMF array access) to the 2.8/im pitch in the 
evanescent coupling region. The coupling region is 
350/im long, before bending gradually to an output pitch 
of 125/im which was accessed by 250/im pitched arrays of 
PMF; the output measurements were alternated between 
access to all odd labeled waveguides and all even labeled 
waveguides. The optical chips were diced from the wafer. 
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Supplementary Information 

Operator for continuous, constant nearest neigh- 
bour coupling. Evanescently coupled uniform waveg- 
uide arrays are modelled with the coupled oscillator 
Hamiltonian 



H 



(10) 



for waveguide propagation constant (site potential) 13 
and constant coupling coefficient C , where /3, C G 
M. In the basis of single position (or mode) states 
{|1)^ , 11)2 , ... |1)„}, the Hamiltonian is represented by a 
tri-diagonal M x M matrix 



Hi , 3 - 



/3, i^j 



(11) 



Unitary evolution of this Hamiltonian after time T = cz 
is given by J7 = e^^^ , where z is propagation distance 
in the waveguide and c is the effective propagation speed 
which is equivalent to a continuous time quantum walk 
on a one dimensional graph [1 . 

The device reported is a uniform array of M = 21 
waveguides. The propagation length z is such that the 
linear ballistic propagation resulting from launching 
light into the central waveguides jlj does not extend 
beyond ten central waveguides which are accessed in 
the experiment. For computing similarity between 
experiment and theory the matrix A is modelled by 
extracting the central 10 x 10 submatrix of a 21 x 21 
unitary matrix U — e*^^, computed using Eq. 10 



with parameters T — 13.9, (3 — Q and C = 0.15; the 
theoretical model assumes the array to be uniform and 
to have zero loss. Both assumptions can give rise to 
S < 100% in our analysis. 



2-particle correlation functions. Quantum interfer- 
ence between two indistinguishable particles — modelled 
by creation operators and annihilation operators Uj 
acting on labeled modes j in a finite set of modes associ- 
ated to system a — is quantified by the correlation func- 
tion and has been derived for photons (bosons) [1] 
and fermions [3] in multi-port networks. Here we re- 
view derivation and adapt for fractional statistics that 
allow continuous deformation between Bose- Einstein and 
Fermi-Dirac statistics. 

Consider two indistinguishable particles that are gov- 
erned by the relations 



«14 



= 



ajCfc - e ^afeOj 



a, a 



(12) 

(13) 
(14) 



for Bose-Einstein (0 — 0) and Fermi-Dirac (0 — tt) statis- 



tics, governing exchange according to the permutation 
group [4]. We can also define particles that undergo a 
fractional phase change Q < <j} < tt with j < k, on ex- 
changing two such particles according to braiding oper- 
ations defined as 'clockwise'. For j > k, anti-clockwise 
braids are defined with — tt < < 0. 

A mode transformation A on a state of two indistin- 
guishable particles created in modes j and k, maps the 
initial state aja| |0) to a linear superposition of two par- 
ticle states according to 



44 



|0) 



A \ - 



AsjAt^kaUl |0) 



(15) 



where A is expressed as a matrix in the mode basis. Cor- 
related detection at r and q will occur as a result of two 
possible events: (i.) the particle at j is mapped to posi- 
tion r and the particle at k is mapped to q with ampli- 
tude ArjAq^k', (ii.) the particle at j is mapped to position 
q and the particle at k is mapped to r with amplitude 
AqjAr^k- Labelling modes such that j < k and r < q 
implies (ii.) is equivalent to at least one exchange of the 
position of the two particles. Here we define (i.) corre- 
sponds to no exchange while (ii.) corresponds to a single 
swap or braid operation according to ( 12p4 1. 



Since the two particles are indistinguishable, quantum 
mechanics dictates outcomes (i.) and (ii.) are indistin- 
guishable and hence the complex probability amplitudes 
Ar.jAq^k and AqjAr,k interfere in the correlation across r 
and q. The correlated detection outcome Tf ^ is therefore 
computed by projecting the state (0| a^ag on the output 



state (15), using the inner product 



arUqa^^al |0) = Sq^^Sr^k + e"''6rjSq^k, 



(16) 



which follows from the relations ( 12p4 ). Combining ex- 
pressions (15) and (16) therefore yields 



" r,q 



qj-^r^k I 



(17) 



for exchange statistics parameterised by cf). Note that 
for fermions (</> = tt) the diagonal elements of the matrix 
T'!^^ — agree with the Pauli-exclusion principle. 



iV-particle correlation functions. Consider N parti- 
cles with exchange parameter (f) and governed by the re- 
lations ([I2 14 ) which are launched into inputs labelled by 
, i'n} of mode transformation A. The 



tuple u = {vi, 1/2, ■■ 
A'^-fold correlation function for detection across output 
A^-tuple fl — {/xi,/X2, ...,^Ar} is derived for multi-particle 
quantum walks with both Fermi-Dirac and Bose-Einstein 
statistics [5]. 

The input state is mapped by arbitrary ^ to a lin- 
ear superposition of A''-particle states at the output of A 
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according to 



N 



N 



n<io)-n 



i=i 




(18) 



Defining i^i <h'2 < ■■■ <i^n and /ii < /i2 < ■•• < Mat we can 
consider the transformation A maps the N particles at 
a to the output fl in A^! indistinguishable possibilities. 
One of these possibilities maps vj to iij for all j which 
we define to not exchange a single pair of particles. The 
remaining N\ — l map Vj to tr^^ for all j for some permu- 
tation (7^ of and are defined to exchange particles with 
the minimum number of pairwise transpositions in clock- 
wise braids that transform to a^. Each braid between 
any two particles makes the combined wave-function ac- 
quire a phase 0; m clockwise braids that are the min- 
imum number required to map v to cFi, therefore cause 
the global wave-function to acquire phase m^. From this 
definition, we obtain an A^-particle generalisation of ex- 



pression ( 16 1: 



Af TV N 

((oi n (n =E n ^^^.-'^.^ (i^) 



where T{a^) is the minimum number of pairwise trans- 
positions (neighbouring swaps) relating fjn to v. Pro- 
jection onto the output state ( [T8| with TV-particles at fl 
uses relations ( 12][T4 ) and ( [T9| for a fractional parameter 
(0 < < tt) to yield 



(20) 



Varying cf) allows continuous deformation between Bose- 
Einstein (</> = 0) and Fermi-Dirac (</) = tt) quantum 
interference. 

Generating A^— partite, A^— level entanglement. A 

quantum algorithm for anti-symmetrising arbitrary prod- 
ucts of n wave functions p2, Pn)^ is presented [6] 
that would use three registers of n, to— digit quantum 
words or "qu- words" (strings of n log2 to qubits) . This 
was then simplified l7j to use two registers of qubits A 
and B] the first encodes the wavefunctions |pi, p2, Pn)- 
The second register is mapped to a symmetrized state 
with 0(n(lnTO)^) local operations followed by 0(n^ In to) 
controlled operations leaving the second register in the 
state 



1 



(21) 



u£Sn 



where a{\,..,N) is a permutation <t G of the 

A^-tuple {1,2, The state Eq. I2l] is then anti- 



symmetrised by a sorting algorithm using O(nlnn) op- 



erations, applying the same swaps (via controlled oper- 
ations) and phase shifts to the first register to antisym- 
metrize \pi, p2, Pn) a- For our purpose, we require only 
the one register B, since the state we wish to generate is 
known and of the form 



= 7^ E e-(-'^)*|a.(l,2,...,iV))(22) 

where (t(1,..,A^) can be decomposed into T(cr) adjacent 
transpositions. 

The approach shown in Fig. [s] represents each i*'' qu- 
word (denoted B[i]) in the register as an TV— level qu- 
dit, for the examples iV = 2, 3, and generates the state 
\'4'N (0)) in an iterative process, using X^ili^ i < 

controlled operations. Here each controlled opera- 
tion is counted as one swap operation of two modes in 
a given qudit, conditional on an excitation in one of the 
N levels of the control qudit. Realisation with linear op- 
tics for example, could use a heralded KLM CNOT gate 
[SlIS], each requiring two extra (ancilla) photons, and an 
interferometric structure of seven beamsplitters. 

As with Refs. IBJ El the N qudits are initialised in su- 
perposition using local operations (mode splitters; M S) 



\^IB[N] ^ 

MS 1 
!• ■ 



lB\l] 



El'=)s[A'] ® E l^)s[iV-l] 
Kk=l ) \k=\ J 



(23) 

B[l] 



Note that this step puts the state in the correct number 
of N\ terms required for the eventual entangled state. 
Using linear optics, this would use one photon for each 



qudit and = 0{N'^) two-mode beam spHtters 

(with appropriate reflectivities). 

The next step in generating \'4>n{4>)) is to introduce the 
phase relationship </> according to 




j{k-i)4, 



\k) 



B[N] _ 



k=l 



(24) 



. fc=i 



which uses ~ 1) = 0{N^) local phase shift oper- 

ations. A set of X^ila^ *X)j=ii = 0{N^) controlled op- 
erations are then applied iteratively to generate |'0Ar (</>)); 
this process can be thought of as performing a series of 
controlled operations on the intermediate states 



1 



V<1~ 
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for each step of the iteration q = 2,...,N — 1. This is 
verified for g = 2 (see Fig. [sjja) and Fig. [5][b)). Next we 
assume inductively the process works for some q ^ r to 
generate lipg^cj))). For q = r + 1, the process would begin 
with the state 



A series of controlled qudit shift operations are then ap- 
plied, such that an excitation of mode j of the i?[r + 1] 
qudit shifts each of the labeled modes {j,...,r} to {j + 
1, ...,r + 1} in the target qudits labelled {-Bi?"], . • .i?[l]}, 
such that the overall state evolves to 



= Ee''^'-'^1j)B[.+il®IV'.(</'))B[.,..,i] 



r' ^ 

■cres. 



(26) 



k(lv,r))j 



__L_ ^ e--^-)^! |a(2, r + 1)) s[r,....i] + e"^ \2) s[r+i] 3, 4, r + l))5[.,...,i] 

■ • • + e'^'-'^^ \j)B[r+i] 2, 3, J - 1, J + 1, r + 1)>b[.,...,i] + ■ • ■ 



(27) 



The tuple {j,l,2,...,j - l,i + l,...,r + 1} is (j - 1) 
adjacent transpositions from the tuple {l,2,...,r + 1}. 
Since the permutation cr(l,2, j — l,j + 1, ...,r + 1) is 
x(ct) adjacent transpositions away from the ordered set 
{1, 2, J — + 1, r + 1}, it follows that the tuple 
{jX1,2, + ...,r+l)} is (j-l+x(cr)) adjacent 

transpositions different from {l,2,..,r + 1}. Hence the 



(a) 



(b) 



Eq. 27 is equivalent to the state 



(28) 

as required. Therefore, by induction we have an iterative 
process that generates \ipNi4')) for arbitrary size q = N, 
requiring O(iV^) local operations, and 0{N'^) controlled 
swap operations. 



1 X 2 mode 
y splitter 





— 






1 


ii— . 



Relative phase 
/ shift d) 



■ Controlled 
mode swap 



Ml] 



<t> 2<t> 
/ 



I A J I 



1x2 7 



FIG. 5: Qudit circuits for geneating \ipi^{(f>)) for (a) N = 2 and 
(b) N = 3. Qudits are initialised in the state |l)g[jv-] ... | ' 

where B[i] denotes the qudit in the register. The green squares 
represent 1 X A; mode sphtters that can be decomposed into k — 1 
two-mode splitters. Blue rectangles represent a relative phase shift 
with labelled phase j<p. Each controlled swap between two modes 
denoted 'x' are conditional on an excitation in the corresponding 
mode denoted denoted by a 
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